Equivalence between fractional exclusion statistics and Fermi liq- 
uid theory in interacting particle systems 

D. V. Anghel^ and G. A. Nemnes^'^ 

^ Horia Hulubei National Institute for Physics and Nuclear Engineering, P.O. Box MG-6, 017126 Mdgurele, Ilfov, 
Romania 

^ University of Bucharest, Faculty of Physics, "Materials and Devices for Electronics and Optoelectronics" Research 
Center, P.O. Box MG-11, 077125 Mdgurele-Ilfov, Romania 



PACS 05 . 30 . -d - Quantum statistical mechanics 

PACS 05 . 30 . Ch - Quantum ensemble theory 

PACS 05.30.Pr - Fractional statistics systems (anyons, etc.) 

Abstract - We develop a method based on fractional exclusion statistics (FES) to describe 
s-dimensional systems of interacting particles. The particles - which are described in the quasi- 
classical approximation - are in an external potential and experience a generic particle-particle 
interaction. We define the FES quasiparticle energies, we calculate the FES parameters of the sys- 
tem and we deduce the equations for the equilibrium particle populations. The FES gas is "ideal", 
in the sense that the quasiparticle energies do not depend on the other quasiparticle levels popu- 
lations and the sum of the quasiparticle energies is equal to the total energy of the system. We 
prove that the FES formalism is equivalent to the standard Landau's Fermi liquid theory (FLT) 
approach and the FES quasiparticle populations may be calculated from the FLT populations by 
making the correspondence between the FES and the FLT quasiparticle energies. The FES pro- 
vides a natural ideal gas description of the interacting particle gas in a generic external potential 
and in any number of dimensions. 



Introduction. — The concept of fractional exclusion 
statistics (FES) was introduced by Haldane in Ref. [1] and 
the statistical mechanics of FES systems was formulated 
by several authors, employing different methods [2H2]. 

The FES was applied to quasiparticle excitations at the 
lowest Landau level in the fractional quantum Hall effect, 
spinon excitations in a spin-^ quantum antiferomagnet 
[TJI31IH], Bose and Fermi systems described in the thermo- 
dynamic Bethe ansatz [SlIMTTj. excitations [1] or motifs 
of spins |121ll3j in spin chains, elementary volumes ob- 
tained by coarse-graining in the phase-space of a system 
[5 |fT4HT8] , interacting particles described in the mean- field 
approximation |19H21] . Thermodynamic properties and 
correlations in generic FES systems in different numbers 
of dimensions or interacting with external fields have also 
been calculated p2t23|23tf27p7t28|28tf3T] . InRefs. [T9ll32] 
it was shown that all the systems with the same, constant 
density of states and diagonal, constant FES parameters 
(FES is manifested within each species) have the same en- 
tropy and heat capacity. This property extends the Bose- 
Fermi thermodynamic equivalence in 2D systems discov- 



ered long time ago [33H35] , to the FES systems [36] . 

A stochastic method for the simulation of the time evo- 
lution of FES systems was introduced in Ref. [37| as a gen- 
eralization of a similar method used for Bose and Fermi 
systems |38j whereas the relatively recent experimental 
realization of the Fermi degeneracy in cold atomic gases 
have renewed the interest in the theoretical investigation 
of non-ideal Fermi systems at low temperatures and their 
interpretation as ideal FES systems [27, 28, 39rHl|. 

The general FES formalism was amended to include the 
change of the FES parameters at the change of the parti- 
cle species [T6ll42fl44] . This amendment allows the general 
implementation of FES as a method for the description of 
interacting particle systems as ideal (quasi)particle sys- 
tems [nnuMzi. 

In this paper we apply FES to a system of particles with 
general particle-particle interaction, l/(|r — r'|), in an ar- 
bitrary external potential, yext(r), described in the quasi- 
classical approximation. We define quasiparticle energies 
which determine our FES parameters and using these we 
calculate the FES equilibrium particle distribution. We 
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compare these results with the standard Fermi hquid the- 
ory (FLT) calculations. 

The structure of the paper is as follows. In the next sec- 
tion we introduce our model and calculate the equilibrium 
particle population in the FLT approach. In the following 
section we implement the FES description by introducing 
new quasiparticle energies and the definition of species. 
These species are related by the FES parameters, which 
we calculate. Using the FES parameters and quasiparticle 
energies we write the equations for the equilibrium parti- 
cle populations. By making the correspondence between 
the FLT and the FES quasiparticle energies correspond- 
ing to the same single-particle states, we show that the 
FES equations are satisfied by the FLT populations. This 
proves that the FES formalism is consistent and suitable 
for the description of such interacting particle systems and 
that the FES quasiparticle populations may be calculated 
using the correspondence between the quasiparticle ener- 
gies. 

The interacting particle system may be described as an 
ideal FES gas in a rather general formalism. 

The model and the FLT approach. We analyze 
a quasiclassical system of particles (bosons or fermions) 
with a generic interaction, in an external potential and s- 
dimensional space. We assume that the particle-particle 
interaction. V{\r — r'|), depends only on the distance 
between the particles, whereas the external potential is 
14xt(r). If we denote by e the energies of the non- 
interacting particles, in the continuous limit the total en- 
ergy of the system is 



where 



E 



e + K=xt(r)+ / dV / de'air',e')n{r',e')V{\r-r'\) 
Jn Jo 

(4) 

are the Landau's quasiparticle energies. 

The equilibrium distribution of particles, n(r, e), is de- 
termined by solving self-consistently Eqs. and (|3]). 
Plugging n(r, e) back into E ^ and ln(Z) ([2]), one can 
calculate the total energy of the system, the grandcanon- 
ical potential and from these, the other thermodynamic 
quantities. We observe that the sum of the quasiparticle 
energies is 

E^ = f d^r [ dee^{e)a{r,e)n{r,e) = E 
Jn Jo 

de' cr(r, e) 



+ - [ d'v [ d^r' [ de 
2 Jn Jn Jo Jo 

xa{v\e')n{v,e)n{v',e')V{\v-v'\), 
which is different from the total energy E. 



(5) 



The FES formalism. — In the FES, the quasiparticle 
energies are defined in such a way that the total energy 
of the system is equal to the sum of quasiparticle energies 
[T7ll45[47|l48] . For this purpose we define new quasiparticle 
energies, 

e~r = e + yext(r)+ /rfV0[e + yext(r)-1/oxt(r')] 



/ d'v f de[e + Kxt(r)]cr(r, e)n(r, e) 
Jn Jo 

TP p pOO poo 

+ - / d'r / d'r' / de de'a{r,e) 
2 Jn Jn Jo Jo 

xa(r',e'Hr,eWr',e')mr-r'|), (1) 



where il is the total volume of the system and n{r, e) is 
the particle population of the level e in the elementary 
volume d^r. In an elementary volume, Sflr, centered at r, 
the density of states (DOS) along the e axis is (50rO'(r, e). 
The logarithm of the grandcanonical partition function is 

ln(Z) = T [ d'r r dea{r,e){[lTn^^\r,e)] 
Jn Jo ^ 

X log[l =F n^^^ (r, e)] ± n'-^^ (r, e) log n^^^ (r, e) | 



n 

e-l-Vext(r)-Kxt(r') 

X J deV(r',e')n(r',e')y(|r-r'|06) 



From Eq. ([6|) we obtain the DOS along the e axis. 



a[r,irie)] = cr(r,e) 



der 



de 



(T(r,e) 



1+ / d^r'0[Fext(r)-yext(r')]mr-r'|) 



(7{r' , e)n[r', e{er 



(7) 



-f3{E~fiN), 



(2) 



where the upper and lower signs correspond to fermions 
and bosons, respectively. 

The standard way to obtain the thermodynamics of the 
system is to maximize Z with respect to n'-^^(r, e). In this 
way we obtain the equilibrium populations, 



To apply FES we split the system into imaginary ele- 
mentary volumes, SQ^ (centered at r), and in each such 
volume the quasiparticle energy axis, e, is split into ele- 
mentary intervals, Se (centered at e). Each such (s -I- 1)- 
dimensional elementary volume, SQ^ x represents a 
FES species. 

From Eqs. ([6]) and ^ and applying the procedure of 
Ref. [T71I3S] we obtain the FES parameters 



(r,e) 



1 



(3) 



o^sn^eisn^re 
oisn^eisnr'^' 

o^snrcisn^, i' 



y(|r-r'|)(7[r,e(e)]<517r (8a) 
^(er-.~;,)n|r-r'|)-'^^^^["(''^)]^ 



de 



xcr[r, e{e)\m^5e = ari.r'i'SG{e), (8b) 
e{i^~i'^,)V{\r~r'\)ar{0)Sn^. (8c) 
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The partition function is the same as the one in Eq. ^ . 



with n'-^^ (r, e) 



= n(±) 



r, er(e)], by definition 



Maxi- 



mizing Z with respect to the population along the e axis, 
we get the equations for the equilibrium particle popula- 



tion 



In 



[lTn(±)(r,e~r)] 



= T 



xV{\r- r'l) [a[r', e{~e,)] ln[l T ^^^^r', S,)] 



de 



, d a{r',e') 



ln[lTn(*)(r',e')] 



d^r'0[Kxt(r') 

|cr(r',0) ln[lT^^ 
, 9cr(r',e 



W(r',0)] 




ln[lTn^^Hr',e')] 



Now we can check the equivalence between the FES and 
the FLT descriptions of the system by substituting Eq. ^ 
into Eq. ^ . By doing so we recover the relation between 
e 



Fig. 1: (Color online) Specific FES and FLT quantities for a 
dimensional system with repulsive Coulomb interactions 
{Nr = 20, Ni = 50, Emax = 25): (a) the quasiparticle density 
of states in the FES description; (b) the position dependent en- 



and e^. 



yext(r')] / de'a{r',e') 



xn(r', e')V{\v - r'|) + B[V,^,{v') - 
' de'a{v',e')n{v',e')V{\v-v'\), 



(10) 



which is in accordance with the definitions (|H) and (jH]). 

We illustrate our model and the relation between the 
FES and the FLT descriptions on a one dimensional sys- 
tem of fermions with repulsive Coulomb interactions - 
y(|r — r'l) = l/|r — r'l, in the absence of external fields. 
The length of the system, J7, is discretized as in Rcf. 

into TVr equal elementary segments, 5r^ = 1, where 
^ = 1, . . . , Nr- In each such elementary "volume" the DOS 
is taken to be constant, ^r^tT(r^, e) = (Tq, and on the e axis 
we define equal consecutive segments between and 
Cmax, 5li = fmax/^e, where i = l,...,iVg - wc choose 
Cmax such that (r, emax) ^ 1 for any r. In this way 
we obtain Nj. x Ni species of particles, 5y^ (E)5ii, identified 
also by a double index, (^, i) [49]. To avoid the singularity 
at the origin of the interaction potential, we consider a 
cut-off distance of 6r^/2 - the choice of the cut-off does 
not affect the procedure. 

The total number of particles in the system is = 
NrCpO'o-, where ej? is the Fermi energy in the noninteract- 
ing system. We set the energy scale of the system by fixing 
ep = 1 and fcsT = 1//3 = 1. 

Figure [TJa) shows the quasiparticle density of states in 
the FES description, a^i = (5r^CT(r^, e,) (Eq. [7]). The den- 
sity of Landau's quasiparticle states may be calculated also 



from Eq. ^ and we obtain = 9{e 



^,min 



ergy shift 



, of the Landau quasiparticle energies; (c) the 



! ^,miii 

FES populations; (d) the particle density, scaled with efcq - 
the dashed line represents the uniform distribution of a similar 
non-interacting system. In the (a) and (c) subplots we plot a 
curve for each elementary segment, (Jrj, ^ = 1 . . . , 20. We use 
the same color for symmetric segments, namely red for 1 and 
N-c, green for 2 and A^r — 1, hlue for 3 and — 2, and black 
for the rest. The chemical potential, /x, is marked by vertical 
dashed lines. 



which is different from zero only for e > e 
imum value of (|4]) is 1^.^^^^^ 



^,min 



The min- 



^ e^fe = 0) = 



The position dependent is plotted in Fig. [T] (b). 



Figure [IJc) shows the populations, 



,(+) 



(|ni). To obtain the Landau's populations ([3]), we simply 
shift the quasiparticle energy from ^ to ef (jU and 
ri^+)[rj, ei(ef )] becomes a Fermi distribution in for any 
^ and with the same /i as for the FES distribution. The 
difference between the distributions n^+-'(r^,ef') in differ- 
ent volumes, ^, being the range of , as shown in Fig. 

ETb). 

The particle density, = ^i^^in^i^ei, is represented 
in Fig. [TJd). Due to the symmetry of the problem, we 
have pair-wise identical populations, n^^i = njv^-^+i^i and 
Pi ^ P-Wr-C+ii for any ^ and i. The largest deviations in 
the particle density occurs for the extremal species, placed 
at both ends of the ID box as an effect of repulsive inter- 
actions. 

We observe from Eqs. ^ and from Fig. [T] that 

the range of ti is [0, oo) in any elementary volume, ^. 
This property of the FES representation allows for par- 
ticle transport along the system at any energy, e or e. 
On the other hand, the range of e^^ is [e^i„in'°°)' with 

mill > for any ^. The presence of creates a 

non-constant energy gap (Fig. [TJ)) for the quasi-particle 
populations and the transport of low energy quasiparticles 
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from one end of the system to the other may occur only 
through tunnehng or higher order processes which involve 
rearrangements of quasiparticles above the energy gap. In 
the FES description these rearrangements are taken care 
of implicitly by the formalism. 

Conclusions. — Wc have formulated an approach by 
which a system of particles with general particle-particle 
interaction, V{\r — r'|), in an s-dimensional space and ex- 
ternal potential Vext(r) is described as an ideal gas of 
fractional exclusion statistics (FES). We have given the 
equations for the calculation of the FES parameters and 
equilibrium populations. 

The FES approach have been compared with the Lan- 
dau's Fermi liquid formalism and we have shown that al- 
though there are differences in the definitions of certain 
quantities like the quasiparticle energies, the physical re- 
sults are the same. The main difference between the two 
formalisms is that in the FES approach the quasiparticle 
energies are independent of the populations of other quasi- 
particle states and therefore the FES gas is "ideal" , with 
the total energy of the gas being equal to the sum of the 
quasiparticle energies. 

We have exemplified our procedure on a one- 
dimensional system of fermions with repulsive Coulomb 
interaction for which we calculated the main microscopic 
parameters, like the quasiparticle energies, quasiparticle 
density of states and energy levels populations. For each 
of these quantities we discussed the similarities and differ- 
ences between the FES and the FLT approaches. 

One practical consequence that appear from our calcula- 
tions is that the solution of the FES integral equations may 
eventually be calculated easier by solving self-consistently 
the FLT equations for population and quasiparticle ener- 
gies, dS]) and g]). 

Another consequence is that while in the FLT formula- 
tion the quasiparticle energies may form an energy gap at 
the lowest end of the spectrum due to the particle-particle 
interaction, in the FES description such an energy gap 
does not exist. For this reason the low energy particle 
tunneling in FLT may be described as higher order pro- 
cesses, involving rearrangements of particles on the higher 
energy levels, wheres in the FES these rearrangements are 
taken into account directly by the manifestation of the 
exclusion statistics. 
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